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Summary. — Research in many areas of modern physics and astrophysics such as, 
e.g., indirect searches for dark matter (DM), particle acceleration in SNR shocks, 
and the spectrum and origin of extragalactic 7-ray background, rely heavily on 
studies of cosmic rays (CRs) and associated diffuse emissions. New or improved 
instrumentation to explore these open issues is ready or under development. A 
fleet of ground-based, balloon-borne, and spacecraft instruments measures many 
CR species, 7-rays, radio, and synchrotron emission. Exploiting the data collected 
by the scientific missions to the fullest requires reliable and detailed calculations 
using a numerical model. GALPROP is the current state-of-the-art numerical CR 
propagation code that has become a standard analysis tool in CR and diffuse 7-ray 
research. It uses astrophysical information, nuclear and particle data as input to 
self-consistently predict CRs, 7-rays, synchrotron emission and other observables. 
This paper reviews recent GALPROP developments and results. 

PACS 95.35.+d - Dark matter (stellar, interstellar, galactic, and cosmological) . 
PACS 98 . 38 . - j - Interstellar medium (ISM) and nebulae in Milky Way. 
PACS 98.70.Rz - 7-ray sources; 7-ray bursts. 

PACS 98 . 70 . Sa - Cosmic rays (including sources, origin, acceleration, and interac- 
tions). 



1. The GALPROP code 

The GALPROP project [1, 2] began in late 1996 and has now 15 years of development 
behind it( x ). The code, originally written in fortran90, was made public in 1998. A 
version rewritten in C++ was produced in 2001, and the most recent public version 54 
was recently released [3]. The code is available from the dedicated website( 2 ) where a 
facility for users to run the code via online forms in a web-browser is also provided. 
The key concept underlying the GALPROP code is that various kinds of data, e.g., 



C ) http: / /sciencewatch.com/dr /crf/2009/09octerf/09octerfStronET / 
( 2 ) http://galprop.stanford.edu 
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direct CR measurements including primary and secondary nuclei, electrons and positrons, 
7-rays, synchrotron radiation, and so forth, are all related to the same astrophysical 
components of the Galaxy and hence have to be modeled self-consistently [4]. The goal 
is for GALPROP-based models to be as realistic as possible and to make use of available 
astronomical information, nuclear and particle data, with a minimum of simplifying 
assumptions. A complete description of the rationale and motivation is given in the 
review [5] . A very short summary of GALPROP is provided below; for details the reader 
is referred to the relevant papers [1, 2, 3, 6, 7, 8, 9, 10, 11]. 

The GALPROP code solves the CR transport equation with a given source distri- 
bution and boundary conditions for all CR species [2]. This includes a galactic wind 
(convection), diffusive reacceleration in the ISM, energy losses, nuclear fragmentation, 
radioactive decay, and production of secondary particles and isotopes. The distribu- 
tion of CR sources can be specified as required ([12] shows an example of the source 
distribution derived from the fit to the Fermi-L AT data). The numerical solution of 
the transport equation is based on a Crank-Nicholson implicit second-order scheme [13]. 
The spatial boundary conditions assume free particle escape. For a given halo size the 
diffusion coefficient, as a function of momentum and the reacceleration or convection 
parameters, is determined from secondary /primary ratios. If reacceleration is included, 
the momentum-space diffusion coefficient D pp is related to the spatial coefficient D xx 
(= (3D a p s ) [14], where 5=1/3 for a Kolmogorov spectrum of interstellar turbulence or 
S = 1/2 for a Kraichnan cascade (but can also be arbitrary), p = pc/Ze is the magnetic 
rigidity. Non-linear wave damping [11] can also be included if required. 

Cross-sections are based on the extensive LANL database, nuclear codes, and param- 
eterizations [15]. The most important isotopic production cross-sections are calculated 
using our fits to major production channels [8, 16]. Other cross-sections are computed 
using phcnomenological approximations [17] and/or [18] renormalized to the data where 
they exist. The nuclear reaction network is built using the Nuclear Data Sheets. Pro- 
duction of neutral pions, secondary positrons and electrons is calculated using the for- 
malism [19, 20] as described in [1] with a correction from [21] or using a parameterization 
given in [22] . Antiproton production uses formalism described in [7] . 

The GALPROP code computes a complete network of primary, secondary and ter- 
tiary CR production starting from input source abundances. Starting with the heaviest 
primary nucleus considered (e.g. 64 Ni) the propagation solution is used to compute the 
source term for its spallation products A— 1, A— 2 and so forth, which are then propagated 
in turn, and so on down to protons, secondary electrons and positrons, and antiprotons. 
To account for some special /3~-decay cases (e.g., 10 Be— s- 10 B) the whole loop is repeated 
twice. The inelastically scattered protons and antiprotons are treated as separate com- 
ponents (secondary protons, tertiary antiprotons). GALPROP includes K-capture and 
electron stripping processes as well as knock-on electrons. 

The 7-rays are calculated using the propagated CR distributions, including a contri- 
bution from secondary particles such as positrons and electrons from inelastic processes 
in the ISM that increases the 7-ray flux at MeV energies [23, 24]. The inverse Compton 
(IC) scattering is treated using the appropriate formalism for an anisotropic radiation 
field [6] with the full spatial and angular distribution of the interstellar radiation field 
(ISRF) [25, 23]. Electron brcmsstrahlung cross section is calculated as described in [9]. 
Gas-related 7-ray intensities (7r°-decay, bremsstrahlung) are computed from the emissivi- 
ties as a function of (R, z, E 7 ) using the column densities of H 1 and H 2 for Galactocentric 
annuli based on recent 21-cm and CO survey data with a more accurate assignment of 
the gas to the Galactocentric rings than earlier versions. The synchrotron emission is 
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computed using a parameterization of the Galactic magnetic field. The line-of-sight 
integration of the corresponding emissivities with the distributions of gas, ISRF, and 
magnetic field yields 7-ray and synchrotron sky maps. Spectra of all species on the cho- 
sen grid and the 7-ray and synchrotron sky maps are output in standard astronomical 
formats for direct comparison with data, e.g., FITS( 3 ), HEALPix( 4 ) [26], Fermi-LAT 
MapCube format for use with LAT Science Tools software( 5 ), etc. 

Also included in GALPROP are specialized routines to calculate the propagation of 
DM annihilation or decay products and associated diffuse 7-ray emission and synchrotron 
sky maps. The routines allow the DM profile, branching ratios, and particle spectra to 
be user-defined and calculate the source functions of the products of DM annihilation 
and 7-ray emissivity. The particles are then propagated as separate species with the 
same propagation parameters as other CRs. The sky maps are calculated using the 
line-of-sight integration of the corresponding emissivities. 

Details of the optimization of the code, linking to other codes (e.g., DarkSUSY [27, 28], 
SuperBayeS [29, 30]) and so forth, can be found at the aforementioned website. 

2. Cosmic-ray propagation in the Galaxy 

Modeling CR propagation in the ISM includes the solution of the partial differential 
equation describing the transport with a given source distribution and boundary con- 
ditions for all CR species. The diffusion-convection equation, sometimes incorporating 
diffusive reacceleration in the ISM, is used for the transport process and has proven to 
be remarkably successful despite its relative simplicity. 

Measurements of stable and radioactive secondary CR nuclei yield the basic infor- 
mation necessary to probe large-scale Galactic properties, such as the diffusion coeffi- 
cient and halo size, the Alfven velocity and/or the convection velocity, as well as the 
mechanisms and sites of CR acceleration. Stable secondary CR nuclei (e.g., 5B) can 
be used to determine ratio of halo size to the diffusion coefficient, while the observed 
abundance of radioactive CR isotopes (4 Be, Al, fyCl, 25M 11 ) allows the separate de- 
termination of halo size and diffusion coefficient [2, 31, 32, 33]. However, the inter- 
pretation of the sharp peaks observed in the secondary to primary CR nuclei ratios 
(e.g., 5B/ 6 C, [2iSc+22Ti+23V]/26Fe) at relatively low energies, ~ 1-few GeV/nucleon, is 
model-dependent. 

Closely connected with the CR propagation is the production of the Galactic diffuse 
7-ray emission, which is comprised of three components: 7r°-decay, bremsstrahlung, and 
IC. Since the 7-rays are undcflectcd by magnetic fields and absorption in the ISM is 
negligible [34], they provide the information necessary to directly probe CR spectra 
and intensities in distant locations, see [35] for a review. However, the interpretation 
of such observations is complicated since the observed 7-ray intensities are the line-of- 
sight integral of a sum of the three components of the diffuse Galactic 7-ray emission, 
an isotropic component (often described as "extragalactic" ) , resolved and unresolved 
sources, together with instrumental background (s). Proper modeling of the diffuse 7- 
ray emission, including the disentanglement of the different components, requires well 
developed models for the ISRF and gas densities, together with the CR propagation [9, 



( 3 ) http://fits.gsfc.nasa.gov/ 

( 4 ) http://healpix.jpl.nasa.gov 
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24]. Secondary CR particles and diffuse 7-rays produced in conventional astrophysical 
processes constitute a background for potential exotic signals (e.g., from DM). 

For details of CR production and propagation the reader is referred to a recent re- 
view [5] . A comprehensive summary of the indirect DM searches in CR and 7-rays can 
be found in [36]. 

3. Recent results 

3'1. Diffuse Galactic and isotropic j-ray emission. - The puzzling "GeV excess" 
relative to the predictions of diffuse 7-ray emission models based on locally measured 
CR spectra [9, 37] was an anomalous signal observed in EGRET data above ~1 GeV. 
It was proposed that the GeV excess results from annihilating DM [38]. This received 
much attention, but a number of conventional explanations were also considered such as, 
e.g., variations in the CR spectra [4, 24]. Paper [39] discusses the sources of systematic 
uncertainties in the EGRET calibration, data handling, and in models of the diffuse 
emission. 

Testing the origin of the GeV excess was one of the early studies of the diffuse 7- 
ray emission by the Fermi-LAT team [40]. The data at intermediate Galactic latitudes 
(10° < |6| < 20°) were used in the study because the diffuse 7-ray emission over this 
region of the sky comes predominantly from relatively nearby CR nuclei interactions 
with interstellar gas. The Fermi-LAT spectrum is well reproduced by the model based 
on local CR measurements and inconsistent with the EGRET GeV excess. Although 
the Fermi-LAT spectral shape is consistent with the model, the overall emission in the 
model predictions using GALPROP was systematically low by 10-20%. This calculation 
employed an a priori model of the diffuse emission, the "conventional" model [9, 24], that 
is based on local CR measurements taken before the Fermi-LAT launch. More detailed 
studies of molecular clouds in the 2nd and 3rd Galactic quadrants [12, 41] show that the 
CR proton spectrum does not fluctuate significantly over a large Galactic volume, which 
supports the reasoning to use the conventional model based on local CR measurements. 
A comparison between the models and the Fermi-LAT and INTEGRAL data in the inner 
Galaxy is discussed in [42]. 

The diffuse Galactic emission presents a strong foreground signal to the much fainter 
diffuse extragalactic emission, which is often referred to as the extragalactic 7-ray back- 
ground (EGB) and generally assumed to have an isotropic or nearly isotropic distribution 
on the sky. The EGB is composed of contributions from unresolved extragalactic sources 
as well as truly diffuse emission processes, such as possible signatures of large-scale struc- 
ture formation, the annihilation or decay of DM, and many other processes [43] . 

The Fermi-LAT measurement of the spectrum of isotropic diffuse 7-ray emission from 
200 MeV to 100 GeV is described in [44]. The isotropic background was found using a 
simultaneous fit of the diffuse Galactic 7-ray emission as modeled using GALPROP, 
resolved sources from the internal Fermi-LAT 9-month source list (using the individual 
localizations but leaving the fluxes in each energy bin to be separately fitted for each 
source), and a model for the solar IC 7-ray emission [45, 46, 47]. The derived EGB 
spectrum is a featureless power law with index 2.41 ±0.05 and intensity /(> 100 MeV) = 
(1.03±0.17)xl0~ 5 cm~ 2 s _1 sr _1 , significantly softer than the one obtained from EGRET 
observations [48] . Note that below 2 GeV the Fermi-LAT spectrum is in agreement with 
the spectrum found from the reanalysis of the EGRET data [49] which was also based 
on GALPROP. Using the fermi-LAT-derived EGB, it was possible to set upper limits 
on the 7-ray flux from cosmological annihilation of DM [50] . 
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Fig. 1. - The luminosity budget of the Milky Way galaxy calculated for a model with 4 kpc 
halo [63]. The percentage figures are shown with respect to the total injected luminosity in CRs. 
The percentages in brackets show the values relative to the luminosity of their respective lepton 
populations (primary electrons, secondary electrons/positrons). 



3'2. Global CR-related luminosity of the Milky Way. - Observations of the diffuse 
7-ray emission from normal galaxies (LMC, SMC, M 31) and the starburst galaxies 
(M 82, NGC 253) by the Fermi-LAT [51, 52, 53, 54] and by the atmospheric Cherenkov 
telescopes [55, 56] show that CRs is a widespread phenomenon associated with the process 
of star formation. The Milky Way is the best-studied non-AGN dominated star-forming 
galaxy, and the only galaxy that direct measurements of CR intensities and spectra are 
available. However, because of our position inside, the derivation of global properties 
is not straightforward and requires detailed models of the spatial distribution of the 
emission. Meanwhile, understanding the global energy budget of processes related to 
the injection and propagation of CRs, and how the energy is distributed across the 
electromagnetic spectrum, is essential to interpret the radio/far-infrared relation [57, 58], 
galactic calorimetry [59], and predictions of extragalactic backgrounds [60, 61, 62], and 
for many other studies. 

Such calculations were carried out in [63]. The luminosity spectra were calculated 
for representative Galactic propagation models that are consistent with CR, radio, and 
7-ray data. Figure 1 shows the detailed energy budget for a model corresponding to the 
middle range of the plausible models. About 1.8% of the total CR luminosity goes into the 
primary and secondary electrons and positrons, however, the IC scattering contributes 
half of the total 7-ray luminosity with the 7r°-decay contributing another half. The 
relationship between far-infrared and radio luminosity appears to be consistent with that 
found for galaxies in general. The Galaxy is found to be nearly a CR electron calorimeter, 
but only if 7-ray emitting processes are taken into account. The synchrotron emission 
alone accounts for only one third of the total electron energy losses with ^10-20% of the 
total synchrotron emission from secondary CR electrons and positrons. 
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Fig. 2. - ID marginalized posterior probability distribution function (PDF) normalized to the 
peak for the diffusion model parameters, with uniform priors assumed over the parameter 
ranges [64]. The cross represents the best fit, the vertical thin line the posterior mean, and 
the horizontal bar the 68% and 95% error ranges, respectively. 



3'3. Constraints on CR propagation models from a global Bayesian analysis. - The 
fully Bayesian approach to the problem of deriving constraints for CR propagation models 
parameters allows one to carry out a global statistical analysis of the whole parameter 
space, rather than be limited to scanning a reduced number of dimensions at the time. 
This is important in order to be able to fit simultaneously all relevant CR parameters 
and to explore degeneracies. While very detailed numerical models of CR propagation 
exist, a quantitative statistical analysis of such models has been so far hampered by 
the large computational effort that those models require. Although statistical analyses 
have been carried out before using semi-analytical models, the evaluation of the results 
obtained from such models is difficult, as they necessarily suffer from many simplifying 
assumptions. 

A full Bayesian parameter estimation has been recently shown to work with a nu- 
merical CR propagation model [64]. Despite the heavy computational demands of a 
numerical propagation code, such as GALPROP, a full Bayesian analysis is possible 
using nested sampling and Markov Chain Monte Carlo methods (implemented in the 
SuperBayeS code [29, 30]). A remarkable agreement was found between the "by-eye" 
fitting in the past [2, 7, 14, 65] and the parameter constraints from the refined Bayesian 
inference analysis (Figure 2) [64]. The posterior mean values of the diffusion coefficient 
D Q = (8.32 ± 4.46) x 10 28 cm 2 s" 1 at 4 GV and the Alfven speed v A \{ = 38.4±2.4 km s" 1 
are in fair agreement with earlier estimates of 5.73 x 10 28 cm 2 s" 1 and 36 km s _1 [41], 
respectively. The posterior mean halo size is 5.4 ± 1.4 kpc, also in agreement with our 
earlier estimated range z^ =4 — 6 kpc [65], although our best-fit value of z^ = 3.9 
kpc is somewhat lower, due to the degeneracy between Dq and z^. However, the well- 
defined posterior intervals produced in that study are significantly more valuable than 
just the best-fit values themselves as they provide an estimate of associated theoretical 
uncertainties and may point to a potential inconsistency between different types of data. 
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